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1. Introduction: AGK cutting rules 

The goal of this paper is to develop a technique that will allow us to calculate more exclusive processes than the total 
cross section in the framework of high density QCD |y, 0, 0]. We will work in the mean field approximation in which 
we have Balitsky-KovchegovM 0] non-linear evolution equation for the elastic high energy amplitude. Except the 
obvious motivation for expanding the calculative power of high density QCD we have a more specific reason for doing 
this paper: to give more transparent derivation of the Kovchegov-Levin equation ||g] for the diffractive dissociation 
cross section in the mean field approximation^. 

The tool which we are going to use, is the AGK gutting rules H. In high energy phenomenology these rules 
are very useful in calculation of the processes of different multiplicity from the expression for the total cross section. 
Since we have a very powerful framework for high density phenomena in QCD, namely, the dipole approach |M|, 
we wish to expand this approach to multiparticle production. The AGK cutting rules, being proved in QCD, will 
allow us to approach such exclusive processes as diffractive production and different correlations in multiparticle 
production processes. For a long time the situation with the proof of the AGK cutting rules has been uncertain (see 
Refs.y, 0, no, O, 12|). At first sight in the leading \og{l/xBjorken) approximation |13| of perturbative QCD the 
scattering amplitude, as a function of particle masses, decreases enough to apply the original arguments of the AGK 
paper M], on the other hand, the main question whether the total cross section and the multiparticle production can 
be described by the same set of diagrams, remains unanswered. The situation became even worse when the explicit 
violation of the AGK cutting rules were found in Refs. [|l^, |ll| |2|. Fortunately, we believe that the mess with AGK 
cutting rules in QCD has been resolved in Rcf. ||lj]. However, before describing the main results of this paper which 
we will use below, we give an introduction to AGK cutting rules explaining what they claim. 



^One of us (E.L.), has heard a lot of complaints from experts that Ref.lql is impossible to understand. 



In simple language, the AGK cutting rules give us the relation between the total cross section at high energy 
and the processes of multiparticle production. In QCD the main idea stems from the unitarity constraint for the 
BFKL Pomeron ||l^, which describes the high energy scattering amplitude in the Leading Log (1/x) Approximation 
of perturbative QCD. The unitarity reads as 



2N{Y-x,v) = \N{Y-x,y)\^ + G,n{Y-x,v) 



(1.1) 



where Y — \og{\/ xsjorken) and [x, y) are the coordinates of the incoming dipole; N {Y; x, y) is the imaginary part of 
the elastic amplitude and the first term describes the elastic scattering (assuming that the real part of the amplitude 
is small at high energy), while the second term stands for the contribution of all inelastic processes. In the leading log 



{l/x) approximation the elastic contribution can be neglected and for the BFKL Pomeron Eq. (1.1) can be reduced 
to the form (see Fig. m: 
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Figure 1: 

ladder. 



(F; X, y) by cut Pomeron while ]\[BFKL ^y . ^^ ^-j ^^y^ ^^ called Pomeron or uncut Pomeron. 

The original AGK cutting rules state that if we 
know the contribution to the total cross section of the 
exchange of any number of Pomerons we can calculate 
the processes with different multiplicity. In Fig. || you 
can see the simple triple Pomeron diagram with the 
AGK coefficients (see Fig. 0-a). The coefficients mean 
that you need to multiply by these coefficients the con- 
tribution of triple Pomeron diagram in the total cross 
section to obtain the cross section of the multiparticle 
production. In Ref. it was proved that these AGK 
cutting rules are correct in QCD for the triple BFKL 
Pomeron vertex. 




Cut Pomeron 
The definition of cut Pomeron through the BFKL 




Double cut 



Single cut 



Diffractive cut 



Fig. l-b 

Figure 2: The AGK cutting rules for the total inclusive processes. Cut Pomeron is defined in Fig. nl and in Eq. ( |l.2| ). In 
Fig. one can see decoding of the Pomeron diagrams in terms of the production processes. 



The exact form of the 


AGK rules look as follows. 


Double Cut: 


^ f (fx2K{xio\xo2,Xi2){2iNBFKLil2) + NBFKL{02)f - 2 7V|j.^i (10)} (1.3) 


Single Cut: 


^Jcfx2Kix,o\xo2,X,2){-^iNBFKL{l2) + A^Bf KL (02))" + 4 7V|j.^i (10)} (1.4) 


Diffractive Cut: 


^ J d^X2K{x,o\xo2,Xi2){{NBFKL{U) + A^Bf KL (02))' - N^pj^^ilO)} (1.5) 


where the vertex of decay 


of one dipole Xio to two dipoles with sizes xq2 and Xi2 is equal to 




°'S T^ f \ \ 0!s xIq 

/^ (.Tin .7;m..T;n) — „ „ ll.fi) 



NBFKL^Xik) is the scattering amplitude of the dipole with size Xik off the target in the leading \og{l/ XBjorken) 
approximation of perturbative QCD (the BFKL Pomeron exchange). It is clear from the above equations that they 
have a simple meaning: the dipole with size xio decays in two dipoles with sizes 2:02 and a;i2 in the case when both of 
them interact with the target with multiparticle production ( double cut of Eq. ( |1.3D ); one interacts with the target 
with multiparticle production while the second one interacts elastically (single cut of Eq. ( |l.4| )); and two produced 
dipoles scatter elastically (diffractive cut of Eq. (|]^)). 

It turns out that difficulties with the AGK cutting rules that has been discussed, stem from the fact that in QCD 
the AGK cutting rules do not work for the situation when we an extra gluon is emitted from the vertices (see Fig. y). 
The difference between these two cases are due to the fact that for multiparticle production without measuring the 





Figure 3: Violation of AGK cutting rules for the processes with additional emission of one particle (gluon). Cut Pomeron is 
defined in Fig. m and in Eq. (1.2). 



gluon, is described by the same set of the diagrams as the total cross section, while for the events with measured 
gluon the set of diagrams turns out to be quite different (see Ref . [^ ) . 

In the next section we will derive the generating functional for multiparticle production in the dipole approach. 
Therefore, for the dipole approach to QCD we will repeat the program that has been worked out for the BFKL 
Pomeron Calculus in zero transverse dimension |L5| . In section 3 we derive the non- linear equation for the generating 
functional and for the cross sections of multiparticle production. These equations will lead to a natural generalization 
of the Kovchegov-Levin equation for the diffractive production cross section (see Ref. g). In conclusions we 
summarize our results. 



2. Generating functional for multiparticle production processes 

2.1 Generating functional for total cross section 

In the mean field approximation (MFA) we take into account only one Pomeron to two Pomeron splitting neglecting 
merging of two Pomerons into one Pomeron. In other words, we consider only 'fan' diagrams. The simple process 



for which such approximation can work, is the deep inelastic scattering with a nuclear target (see Ref.lq]). In the 
dipole approach, one Pomeron to two Ponieron splitting reduces to the decay of one dipole to two dipoles with the 



probability of this decay given by Eq. (1.6)|8|. The simplest and the most transparent technique to incorporate this 
decay is the generating functional which allows us to reduce the calculation of the high energy elastic amplitude to 
consideration of a Markov process, The advantage of this approach is the fact that it takes into account both t and 
s-channel unitarity as we will discuss below. 

The generating functional is defined asP, |l6| 

/n 
P„(y-2/;ri,...,r„) JJ uin)d\, (2.1) 

where u{ri) is an arbitrary function of r.^ and bi. Pn is the probability density to find n dipoles with sizes ri, . . . , r„ 
at rapidity Y ~ y. ^ For functional of Eq. (2T) we have two obvious conditions: 



• at Y = y Pn = for n > 1 and Pi = 5{f — fi). In other words, at y = Y we have one dipole of size r 

Zo{Y-y = 0;{u}) = u(r) ; (2.2) 

• at u = 1 

Zo{Y-y;{u})\u=i = 1. (2.3) 

Eq. (^.3|) follows from the physical meaning of P„ and represents the conservation of the total probability. 



For probabilities P„ we can write the following equation 

_ dPnjy -y-,ri,...,r,,...,rn) ^ 
dy 

- n ^ 

- ? E / '^'^' ^ l*^^; ^'' I^^J- - '^^l) ^" iY-y;n,...,rj,..., r„) 

n-l 



(2.4) 



— ^ / d^r'd^fj5{fj - Ti + f.n)5{r' - f„) K (fj;r', \fj - r'| j _P„_i (Y - y;ri, . . . ,rj, . . . , r„„i) 



Eq. (2^) describes a typical Markov process: two terms of this equation has simple meaning of increase of the 



probability to find n-dipoles due to decay of one dipole to two dipoles (birth terms, the second term in Eq. (2.4)) and 



of decrease of the probability since one of n-dipoles can decay (death term, the first tern in Eq. (2.4)). Multiplying 



by product Jli=i ''^(^i) ^-i^d integrating over r^ we obtain the following linear equation for the generating functional 
^^^l(L_^lMl= If/ d\d\'K{r;r',\f-P\){-u{r) + u{r')u{\P-r\)}^^Zo{Y-y;{u}){2.5) 

Here we use notation S/Su{r) for the functional derivative. Using the generating functional we can calculate the 
scattering amplitude N{Y, r, b) as follows [ [l6| 



^The dipole (xi,yi) with coordinates Xi for qua rk a nd yi for antiquark can be characterized by the dipole size fi = Xi — ffi and 
bi = -kixi + i/i) ■ For simplicity we suppress in Eq. (2.1) and below the coordinate bi. For the scattering with the nuclear target we can 



consider that impact parameters of all dipoles are the same bi = b (see Ref . [pi ) . 
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u— 1— 7 
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where 7(r) = ^(i^o, J', b) is the dipole scattering ampUtude at low energy (Vo). The notation Zq {Y — Yq; {u}) 



(2.6) 

(2.7) 

\u{r) 



means that the initial condition of the generating functional is given by Eq. (2.2). In a more general case we can define 
the scattering amplitude for any arbitrary initial condition F{u{r)), for example, if we choose F{u{r)) = M^(r), we 
start with two dipoles at the initial rapidity Y. In Eq. ( p.6|) we assume that the low energy amplitude of interaction 
of n dipoles with the target 7„ is equal to 7„ = 7", which means that dipoles interact with the target independently 
(without correlations). 



Eq. (2.5) being a linear equation with only first derivative has a general solution of the form Zq (Y — y\ {u\) = 
Z{) {{u{Y — J/)}). Inserting this solution and using the initial condition of Eq. (2.2) we obtain the non-linear equation 



dZ„ (y - y- {u}) r''') _ as 
dy 2-K 



d\'K(r-r', \r- r'\) {-Zo (Y - y,M) r*"-' + Zo {Y - y,{u}) ^^^'^ Zo [Y - y,{u}) p'(l''-'-'l)} 

(2.8) 



It is easy to see that using Eq. (2/7) for the amplitude Eq. (2_^) can be re- written as the Balitsky-Kovchegov equation 
[0, Bi in the form 



dN {Y - y, r, b) 



as 

2tt 



Sr' 



X |iV(y-j/,r',6) + N{Y ~y,\r-?lb) - N{Y -y,r,b) - N{Y -y,r',b) N{Y - y,\f - P\,b)\ 



The evolution equation in its linear form Eq. (2.5) has two advantages in comparison with the non-linear equations 



Eq. (2.6): it has a simple statistical interpretation and can be solved with arbitrary initial condition while non-linear 



equation are correct only with specific initial condition of Eq. (2.2). The statistical interpretation allows us to reduce 
the problem of s-channel unitarity to the conservation of probabilities at each level of rapidity y which is included 



in Markov chain equation (see Eq. {2A) ). Since the technique is equivalent to summing of the Pomeron diagrams 
the i-channel unitarity is preserved, but we need to remember that we must sum all Pomeron diagrams to fulfill 
i-channcl unitarity including the so called Pomeron loops. 

The equivalence between the generating functional approach and the Pomeron Calculus becomes clear if we 
introduce a new functional 

N{Y-y,{j}) = 1- Zo{Y-y;{u})\u=i-j (2.10) 



for which we have the equation 

d{Y-y) 2tt J V ' ' ' 

Starting with the initial condition 



1) {7('-') + lir- r') - 7(r) - i{r') 7(r - ?)] j^N (J - y, {7}) 



(2.11) 



^(i"o;{7})=7(0 



(2.12) 



found from Eq. (|2.2|) and Eq. (p.7|), one can easily built the Pomeron Calculus iterating Eq. (2.11 



2.2 Generating functional for multiparticle production: general ideas 



In the spirit of Ref. [h5| we wish to introduce a generahzation of Eq. (2.1), namely, 



oo „ 71 m 

Z{Y-y;{u},{v}) = ^ P;^ {y;ri, ■ ■ ■ ,rn;n, . . . ,r„,) H u{n)l[ v{rk)d^nd^n (2.13) 

n=l,m=Q •' i=l k=l 

where P™ is the probability to find n uncut Pomerons and m cut Pomerons. In our discussion we use the fact 
that colorless dipoles are proper degrees of freedom for partonic description of the BFKL Pomeron. To treat this 
probability in terms of the dipole approach we need to introduce two typical moments of time for the processes of 
multiparticle production: r = when the interaction with the target occurs, and r = oo where our detectors for 
particles are located. Having in mind these two moments of time we can treat -P™ as the probability to have at 
rapidity y n dipoles with sizes ri, . . . , r„ at t = 0, which do not survive until t = oo and cannot be measured, while 
the dipoles with sizes ri, . . . , r™ reach r = oo and can be caught by detectors. 

The first boundary condition for the generating functional is obvious 

Z{Y-y;{u},{v})\u=i,v=i = 1 (2.14) 

and it follows from the physical meaning of P™ and represents the conservation of the total probability at any 
rapidity. 

Generalizing Eq. (vM) we can introduce unintegrated over impact parameter cross section for multiparticle 
production as 

M{Y-Yo;r,b) = (2.15) 

J^ ( _^\n-{-m-\-l 






'^ — ' ^n ^''. I n A„(,,) >.(..) j-i^ j-i^ U(r-r.);{„))i:<:>,... 



^-^ n\m\ , ^ 

n— l.rn— ^^ i— 1 k—1 



z(r-ro;M,M)i:L1-,,„=i-,.. 



The total cross section is obtained from Eq. (2.15) substituting 27 (r) = "fmir), where 7(7") is the low energy (Yq) 



amplitude of the dipole interaction with the target. 27(r) — ^in{f) follows directly from the unitarity constraint of 



Eq. (1.2) at low energies. We want to remind at this point that we consider Iq is so large that asY^ « 1 and we can 
use the Born approximation of perturbative QCD for estimates of 7^ (r) . Using new notation for elastic amplitude 
A^ (see Eq. ( |2.6| )) and generating functionals Zo and Z (see Eq. (2.1) and Eq. ( 2.13| )), we can re- write Eq. (1.1) in 



the form 

2iV(r-2/,r,6) = 2(l-Zo(r-y;M)i:L1-^) = |7V(y-y,r,6)p + G.„ = 1 - Z {Y - y-{u}Av})\lt\-^,,=,.2^ 

(2.16) 
which translates into boundary condition 

Z[Y-y- {u}, {v}) |,=2„-i = 2 Zo (y - y; {u}) - 1 (2.17) 

It is useful for our further discussions to introduce a generating functional for multiparticle production with an 
arbitrary initial condition 

M{Y - y- {7}, {7„}) ^l-Z{Y-y- {u}, {v})U=i-^,.=i^^,^ (2.18) 



for which boundary condition Eq. (2.17) takes form of 



M{Y - y- {7}, {7m})l7>„=2^ = 2iV(r - y; {7}) (2.19) 




Y 0=0 



Figure 4: The examples of Pomeron diagrams tiiat 
contribute to the process of single diffraction in the 
mean field approximation. The Pomeron, crossed by 
the dotted line,is the cut Pomeron. 



The generating functionals Eq. (2.13) and Eq. (2.18) should 
be supplemented by some initial condition for calculating physical 
observables. The initial condition depends on the process we want 
to calculate. As it was already mentioned, the linear equation for 
the generating functional will not depend on the initial condition. 

As a simple example we calculate the cross section of single 
diffraction the initial condition for the generating functional has 
the following form 

Z{Y;{u},{v}) - v{r) (2.20) 

which means that we consider diagrams which start with one cut 



Pomeron (see Fig. ^. 

The cross section of the single diffractive production with mass smaller than \n{A'l'^ /tv?) < Y — Y' where m is 
the nucleon mass, is equal to (see more details in Rcf. |l5| ) 

asAYX) = 1 - Z{Y-Y'Auh{v})\trWNiY'-,rM.,.ir) = l-miY'-,r,b) (2-21) 

where N{Y'; r, b) reflects the fact that each Pomeron that cross line Y' can develop its own tree 'fan' of Pomerons. 

The substitution v{r) = 1 — N'^{Y';r,b) is explained as follows. Any dipole that corresponds to v{r) survives 
until T = oo ( is present on the unitarity cut) and thus can scatter both elastically and inelastically by 2N{Y'; r, b) — 
\N{Y'; r, 5)^ + G™ (see normalization of -fin = 27). However, we are interested in the single diffractive process with 
no particle production (inelastic scattering) below some rapidity Y'. That is the reason why we retain only the elastic 
part of the total cross section of a dipole present on the cut, namely, N{Y'; r, b)^. 



2.3 AGK cutting rules and vertices for dipole decays 



We wish to write the equation of Eq. (2.4)-typc for P™, but first we need to determine the vertices of decays of one 
v-dipole to two u-dipoles; one u-dipole to one w-dipole and one w-dipole; and one w-dipole to two w-dipoles. 

We will use the AGK cutting rules of Eq. (|l.3|)-Eq. (1.5) to find out these vertices. According section 2.1 we 
expect the following rules for the different transition 
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(2.22) 
(2.23) 
(2.24) 
(2.25) 



where p denotes the cut BFKL Pomeron. At first sight, comparing Eqs. (2.22)-(p.25|) with the AGK cutting rules 



of Eqs. (1L3D-(|1^ for P ^ P + P amip ^ P we have 



P^P + p 
P^P 



d r2K{rio\ri2,ro2) { | {ltn{ri2) 



lin{r02)) 



as_ 

2tt 

— I d'^r2K{rio\ri2,ro2){'-finiri2)+lin{ro2) -lin{rio)} 



1 i,,„,| '^MM-rlMf 

2 J S -fin 



SM{y;{j},{-f,n}) 

5lm 



(2.27) 



However, the AGK rules developed in flS] were written for amplitudes, which are the solutions to the linear (BFKL 
and a generalized form of BFKL) equations. In the present study we are interested in the non-linear evolution and 
thus all quadratic terms of the same argument (7j^„(12) etc.) are by definition absorbed in the corresponding linear 
terms, in other words for the generating functional M{y; {7}, {7i„}) functions 7(r) and 7m (r) are arbitrary and any 
quadratic term of the same argument can be absorbed in the definition of the corresponding linear term. This means 
that the proper way to account for the the AGK cutting rules is as follows 



fP^fP + fP 
P -f JP + P 
P -f P + P 

P -> P + P 
P^P 



as_ 

2tt 



(Pr2 K (rio|ri2, ro2) 7m(''i2) 7m(^02) 



^Af(y;{7},{7m}) 



-— / d?r2K{rw\ri2,rm)l{ri2)l{ro2) 



- 7m (^02)' 
(5Af(y;{7},{7m}) 



(2.28) 



-2^ /d2r2 if (rio|ri2,ro2){7» (^12) 7(^12) + 7^(^02)7(^02)} ^-^^^f^^^^2.2%) 
27r J 7i„(rio) 
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27r 



5^in{riQ) 
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27r 



(fr2 K (ri(^\ri2,rQ2) {'yin{ri2) + Imirm) -7m(rio)} 



JM(^;{7},{7m}) 



(fr2K{rio\ri2,ro2) 7(''i2) 7(^02) 



SM{y;{j},H„}) 



Slirio) 



as 

27r 



-— / d^r2-ftr(rio|ri2,ro2){7(^i2) +7('"02) -7(''io)} 



(5Af(y;{7},{7m}) 



(2.30) 

(2.31) 
(2.32) 
(2.33) 



From these equations we can easily to build the Pomeron Calculus for cut and uncut Pomerons using the simple 
values: S/S"f{r) and S/d"fin{r) are the annihilation operator for uncut and cut Pomerons while the multiplication by 



7 and Jin leads to a creation of uncut and cut Pomerons. It should be stressed that Eq. (2.28)-Eq. ( p.33| ) give a 
direct generalization of the equations for total cross section. In the latter we take into account the sum of all cuts. 
This sum results in the equation in which transitions of Eq. (2.32) and Eq. ( ^.33 ) remain, but with opposite signs. 
It is clear that summing all the cuts we obtain Eq. (2.11). Therefore, we check a selfconsistence of our approach this 
way at the end of the present analysis. 



2.4 Linear functional equation 



Eq. (2.28)-Eq. (2.33) allow us to write the linear evolution equation for N of Eq. (2.15). It has the following form 
aM(y,{7},{7.„}) 



dy 



a.s 



2-K 



—— d r2K {rio\ri2,ro2) < {"/(ri2) + jiro2) — "/{rio) — "f (r 12) "f{r 12)) 



SM{y,{j},{y,„}) 
Sj{rio) 



+ 



(2.34) 



+ (7in(ri2) + 7in(?-o2) - Jinirio) + 7in ('"12 ) 7™ (''02) - 2 7i„ (ri2) 7(ro2) - 2 7i„(ro2)7(^i2) + 2 7(ri2)7(r()2)) 



SM{y,{j},{'y,„}) 
Sjinirio) 



The fact that we have different signs in front of terms jf^, 7m7 and 7^ (see for example Eq. ( 2.29 )), makes the 



probabilistic interpretation very questionable. However, it turns out that we can reduce Eq. (2.34) to a very simple 
equation with a very transparent probabilistic interpretation if we go back to the generating functional Z (y; {u}, {v}) 
where m = 1 — 7 and v = 1 — jin- As it was shown in Pq] , Z (y; {u}, {v}) can be written as a functional of two other 
functions such that Z (y; {u}, {v}) — Z {y; {u}, {£,}), where we defined a new function ^(r) = 2u{r) — w(r). In terms 



of functions u and ^ the linear equation Eq. (2.34) has a simple form 



dZ{y-{u}Ai}) 
dy 



as 



— / d r2iir(rio|ri2,ro2) 



(2.35) 



(u(ri2)u(ro2) - w(rio)) r— p r h {i(rxi)i(rQi) - C(rio)) 



Su{rio) 



S^rio) 



Eq. ( 2.35 ) has a transparent probabilistic interpretation, since the second term at the r.h.s. of the equation is of 
the same structure as the first one and describes the Markov chain for the decay of the dipole which corresponds to 
2u — V. This physical meaning becomes clear if we recall that in the generating functional approach the functions u 
or V correspond to the creation operator for dipole that docs not or does survive till time t — oo, respectively (for 
corresponding annihilation operators we have S/6u and 6/Sv). 

The solution to Eq. (|2.35|) with the initial condition of Eq. ( ^2.20 ) is found as 



^(y;M,{a)r"'''^"^*^^ = 2Zo(2/;M)rM - Zo (y; U}) l^^''^ 



(2.36) 



One can easily check that the solution Eq. ( 2.36| ) satisfies boundary conditions Eq. (2.14) and Eq. (2.17). It is 
very instructive to compare this with the explicit form of the solution in the toy model found in Ref. |lq| . 



3. Non-linear evolution equation 



In this section we rewrite the linear functional first order differential equation Eq. (2.35) as the non- linear equation 
for the scattering amplitude using the initial condition for the generating functionals Zq and Z, given by Eq. ( p.2[) 
and Eq. ( 2.20 ), respectively. The easiest way is to go back to Eq. ( 2.34 ) and using the initial condition for N{y; {7}) 



(see Eq. (2.12)) and M{y; {7}, {7m})) to write the non-linear equation in analogy with transition from Eq. ( |2.5| ) to 
Eq. (2.8). The initial condition for the cross section M (y; {7}, {'jin}) is easily obtained from Eq. ( ^.18 ) and Eq. (2.20) 
in the form of 

M{y = 0;{^},Hn})=%n{r) (3.1) 



Plugging Eq. (2.12|) and Eq. (3.1) into the linear equation Eq. (2.34) we get the non-linear equation 



{M(y;{7},{7™})r»"(^-) + M(y; {7}, {7„})r'"('--) - M(y; {7}, {7,„})r»"('^-) 

+M(y; {7}, HnW-^'-'^^Miy; {7}, HnW^"^''^"^ ~ 2M{y; {7}, HnW'"^''''^ N{y; {-/W^''^"^ 

- 2iV(y; {7})r('--)M(y; {7}, HnW^'^'^'"^ + 2iV(y; {jW^'-'^^Niy; MW^"-^"^} 



(3.2) 



With the help of Eq. (^), Eq. ( |2.10| ) Eq. ( |2.15| ) Eq. ( |2.18| ) we recast Eq. (|j) in a famihar form of 
dM{y;rio,b) _ as /" ,2 



dY 



as 

27r 



d r2 K {rio\ri2,ro2) {M{y;ri2,b) + M{y;r2o,b) - M{y;rio,b) 



(3.3) 



-M(y;ri2,5)M(y;r2o,&) - 2 Af (y; ri2, 6)iV(y; r2o, 6) - 2 iV(y; ri2, 5)M(y; r2o, 6) + 2N{y;n2,b)Niy;r20:b)} 



Eq. ( |3.3| ) describes all multiparticle production processes and has the same form as the equation for the cross section 
of diffractive production found in Ref. M . 



The difference between various processes manifests itself only in the different initial conditions. For example for 



the diffractive production the initial condition has the following form (see Fig. W and Eq. (2.21)) 



M {Y ^ Y'; r, b) = N^ (Y' , r, b) (3.4) 



where iVp is the amplitude of elastic scattering at rapidity Y = Y' . Eq. ( p.4[ ) can be directly obtained from 
the corresponding condition for the generating functional Eq. ( ^.21 ) noting that M (Y ~Y'] r, b) has meaning of the 



total cross section. This way we show that we can describe a process either in terms of the linear equation for the 
generating functional or a non-linear evolution equation for the scattering amplitude (total cross section). It should 
be noted that the generating functional approach is more general and allows more freedom in the choice of the initial 
condition. 

If we want to find a cross section with the fc-recoiled nucleons in dipole nucleus interaction we need to calculate 

a^^'HY^r) = ln^,4rO(^^;^^Af(y;{7},{7.n})r-(^))k^^^ (3.5) 

where 7(r) is the low energy elastic amplitude and 7mi('') = 27 (r) at low energy (Iq)- 

This cross section can be found as the solution to non-linear equation with the initial condition 

aW(ro,r;6) = 1 e^^^^^-'-^'') f] (ro,r; 6) (3.6) 

where VL(Yo,r]b) = <^ (\i^o\e--pmio\S^ = Fo,r)TA (&) with o-dipolc-proton(^ = Yo.r) is the cross section of dipole 
- nucleon interaction at low energy being equal to o'j^ipole-Droton(^ ~ ^o,?") = 2 J (Pb' NQ(Yo,r,b') . The function 
Ta (b) is the optical width of nucleus which gives the number of nucleons at given value of impact parameter b. We 



assume that Yq is large enough to use the Glauber approach. The derivation of Eq. (3^) is given, for example, in 
Ref. [§1. 



4. Conclusions 

The main result of the paper is two equations for the generating functional for multiparticle production: the linear 



equation Eq. ( 2.34 ) and Eq. ( 2.35 ), and the non-linear equation Eq. (3^). The linear equations have an advantage 
of being correct for any initial condition, while the non-linear equation describes the process which starts from the 
exchange of a single Pomeron at low energies. 

The non-linear equation has the same form as the equation for diffractive production that has been proved in 
Ref. and confirms in Refs. ||l^, |l^, |l^. We hope that here we give a more transparent and physically motivated 
derivation for diffractive production and generalize the approach to other processes of multiparticle generation. It 
should be stressed that the processes that we considered here are totally inclusive in the sense that we do not 
measure a particular particle in these processes. For example, our equations cannot describe the multiparticle 
inclusive correlations since we do not have AGK cutting rules for vertices with the emission of gluon (see Fig. 0). 

In this paper we consider only the multiparticle processes in MFA, but we hope to use Mueller-Patel-Salam-Iancu 
approach pc|] to calculate these processes taking into account the Pomeron loops in spirit of the approach suggested 
in Refs. |l|g|l. 
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